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The electromagnetic transition between the almost degenerate 5/2 + and 3/2 + states in 229 Th is 
deemed to be very sensitive to potential changes in the fine structure constant a. State of the art 
Hartree-Fock and Hartree-Fock-Bogoliubov calculations are performed to compute the difference 
in Coulomb energies of the two states which determines the amplification of variations in a into 
variations of the transition frequency. The kinetic energies are also calculated which reflect a possible 
variation in the nucleon or quark masses. A generalized Hellmann-Feynman theorem is proved 
including the use of density-matrix functionals. As the two states differ mainly in the orbit occupied 
by the last unpaired neutron the Coulomb energy difference results from a change in the nuclear 
polarization of the proton distribution. This effect turns out to be rather small and to depend on 
the nuclear model, the amplification varies between about —4 x f0 4 and +4 x f0 4 . Therefore much 
more effort must be put into the improvement of the nuclear models before one can draw conclusions 
from a measured drift in the transition frequency on a temporal drift of fundamental constants. All 
calculations published so far do not reach the necessary fidelity. 

PACS numbers: 06.20. Jr,21. 60. Jz, 27.90. +b 



I. INTRODUCTION 

Measurements of a suspected temporal variation of the 
fine structure constant by means of atomic transitions 
have reached an limit for da/a of less than 10 _16 yr _1 
[1-6]. 

A tempting idea to increase the sensitivity limit is to 
use the transition between two states with very differ- 
ent Coulomb energies, because the drift in transition fre- 
quency 
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is given by an amplification factor A times the drift in 
a. Using the Hellmann-Feynman theorem A is the ra- 
tio of the difference AVc m Coulomb energies of the two 
states involved in the transition divided by the transition 
frequency to. A promising candidate for that is the transi- 
tion between the 3/2 + isomeric state of the nucleus 229 Th 
to its 5/2+ ground state [7, 8]. Recent measurements 
yield to = 7.6 ± 0.5 eV [9] which on nuclear energy scales 
is an accidental almost degeneracy. Typical Coulomb en- 
ergies for this nucleus are of the order of Vc « 10 9 eV so 
that even a small difference AVc could result in a large 
amplification. 

A big drawback of this idea is that the Coulomb en- 
ergies cannot be measured but have to be calculated 
with sufficient accuracy in a nuclear model. Differences 
in charge radii and quadrupole moments of the nuclear 



states are in principle experimentally accessible. Their 
measurement would reduce the uncertainty. 

In a simplified picture the two states differ by the oc- 
cupation of the last neutron orbit. The change in the 
Coulomb energy is thus due to a modified neutron dis- 
tribution in the excited state that will polarize via the 
strong interaction the proton distribution in a slightly 
different way than in the ground state. 

In this paper we investigate in how far state of the art 
nuclear models can provide reliable answers. For that in 
section II we revisit first the Hellmann-Feynman theo- 
rem proposed by H. Hcllmann [10] and also by R. Feyn- 
man [11] and show that it is also valid in a wider con- 
text including density-matrix functional theory employed 
in nuclear physics. After introducing the nuclear mod- 
els in Sec. Ill we discuss in Sec. IV the amplification of 
the temporal drift of fundamental constants like the fine 
structure constant or the nucleon mass when monitoring 
transition frequencies. 

For the nuclear candidate 229 Th we perform in Sec. V 
self-consistent nuclear structure calculations and deter- 
mine via the Hellmann-Feynman theorem the derivatives 
of the transition frequency with respect to the fine struc- 
ture constant and the nucleon mass. A critical assess- 
ment is made on the predictive power when calculating 
observables other than the energy. Finally we summa- 
rize. 
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II. HELLMANN-FEYNMAN THEOREM 
REVISITED 

In this section we show that the Hellmann-Fcynman 
theorem holds not only for eigenstates of the Hamilto- 
nian but also for all stationary solutions in approximate 
schemes, provided they are variational. This is of impor- 
tance as the models we are using are based on density- 
matrix functionals. 

Let the Hamiltonian H(c) depend on a set of external 
parameters c = {ci, C2, . . . }, like the strength of an in- 
teraction or the mass of the particles. The steady state 
solutions of the time-dependent Schrodinger equation are 
given by the eigenvalue problem 



H(c 



(c)| tf n ;c) = E n (c)\ *„;c> 



(2) 



Both, the energies E n (c) and the eigenstates | c ) de- 
pend on c. When discussing intcrmolccular forces in 1933 
H. Hellmann [10] and later in 1939 R. P. Feynman [11] 
showed that a small variation of an external parameter 
in the Hamiltonian (distance between nuclei) leads to a 
change in the energy given by: 
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The proof hinges on | c ) being an eigenstate of the 
hcrmitian H(c): 
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Inserting (2) in the terms enclosed in square brackets 
leads to a cancellation of these terms and we are left with 
Eq. (3), the Hellmann-Fcynman theorem. The partial 
derivative of any eigenvalue E n {c) with respect to a pa- 
rameter Ci is equal to the expectation value of the partial 
derivative of the Hamiltonian calculated with the corre- 
sponding eigenstate | \I>„;c). 

We now show that the same statement holds for any 
extremal point in an energy functional. Let £(c,x) be 
the energy of a physical system that depends on external 
parameters c and on a set of variational parameters x = 
{x\, X2, ■ ■ ■ } which characterize the state of the system, 
x may also represent a set of functions in which case 
partial derivatives are replaced by functional derivatives. 

For example, in the Hartrec-Fock approximation x 
would be the set of occupied single-particle states that 



form a Slater determinant. In an energy density func- 
tional x could be the local density p(r), and so on. 

Steady state solutions x' n )(c), n = 0, 1, 2, . . . of the 
system are obtained by the condition 







7T- c ' x 



(5) 



At the stationary points the energy assumes the values 

£„(c) = 5(c, X W(c)), n = 0,l,2,.... (6) 

Both, the energies and the parameters x( n ) (c) character- 
izing the stationary states depend on the constants c. In 
the ground state given by x(°) the energy £(c,x) is in an 
absolute minimum with respect to variations in x, while 
the other possible solutions x",n 7^ 0, represent saddle 
points. 

A variation of the external parameters at the station- 
ary points leads to 
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Due to the stationarity condition (5) the second part in 
the square brackets vanishes so that we obtain for sta- 
tionary solutions 



£M»>(c)) • 



(8) 



The derivative of the energy at the stationary solutions 
is just the partial derivative of the energy functional with 
respect to the external parameter calculated with the sta- 
tionary state. This generalizes the Hellmann-Fcynman 
theorem (3). 
Taking 
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and I x) = y^x k \ k) , (9) 



where | k ^ denotes some fixed basis and doing a variation 
with respect to the parameters x = {x\,X2, ■ ■ ■ } yields 
the time-independent Schrodinger equation (2) and the 
original Hcllmann-Feynman theorem (3) as a special case 
of (8). _ 

It is interesting to note that the Hellmann-Feynman 
theorem even applies to every solution |\I/(i)) of the 
time-dependent Schrodinger equation 



i-\nt)) = H{c{t)) |*(t)) 



(10) 
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For the total time derivative of the mean energy, which 
is the expectation value of the Hamiltonian, one obtains 

dE _ d (*(*) \H(c(t)) |*(t)) 
~dt ~ 



dt 
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Here we used the fact that a hermitian Hamiltonian does 
not change the norm of the state |\l/(f)). Inserting 
the time-dependent Schrodinger equation (10) into the 
expression in the square brackets lets this term of the 
equation vanish so that the time-dependent Hellmann- 
Fcynman theorem reads: 



dE 
~dt 
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dt 



(12) 



Again the variation of the energy is given by the expec- 
tation value of the derivatives of the Hamiltonian with 
respect to the external parameters. See also Ref. [12]. 

The generalized theorem (8) is even valid for classical 
mechanics, just take x to be the canonical variables p 
and q. The equations of motion are 

dp dE dq d£ 

- = --(c,p,q) and - = -(c,p,q) . (13) 

At the stationary points the partial derivatives of the 
Hamilton function £(c,p, q) vanish and one obtains again 
Eq. (8). 



III. MODELS 

In this section we discuss briefly the Hartee-Fock (HF) 
method when a Hamiltonian is used, the extension to HF 
with density-matrix functionals and the inclusion of pair- 
ing correlations. We show that in all cases the general- 
ized Hellmann-Feynman relation holds. We also explain 
in short the quantities discussed in the section containing 
calculations for 229 Th. 



A. Hartree-Fock with Hamiltonian 

In the HF approximation one uses a single Slater de- 
terminant 



t J 



l 2 ■ 



(14) 



as the many-body trial state. The creation operators aj, 
that create the occupied single-particle states | </>„ ^ , 



4|0> = |^) 



£4 



(15) 



are represented in a working basis | i \ = cj | ). Thus, 
the expansion coefficients D iu represent the set x = 
{Di v ; v = 1, . . . , A, i = 1,2,...} of variational parame- 
ters. 

The general definition of the one-body density operator 



^2 | % > pik < k 



(16) 



i . A- 



is given in terms of the expectation values of c k Ci, where 
the creation operators c- create the single-particle basis 
| % ). In the HF case: 



Pik 
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The energy of the HF Slater determinant can be ex- 
pressed in terms of the idempotent (p 2 = p) one-body 
density as 



£hf[c, p 



( 1>HF | H{c) | ^HF ) 
( *HF | *HF ) 

C ^ Pfi + 4 E V ihjl( c ) PjiPlk 
ijkl 

Vikm,jln( c ) PjiPlkPn 
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(18) 



where i,j (c) denotes the matrix elements of the kinetic 
energy, u^/^c) the antisymmetrized matrix elements of 
the two-body interaction, Vik m ,jln(c) of the three-body 
interaction, and so on. The dependence on the pa- 
rameter set c that includes the nucleon masses, cou- 
pling strengths, interaction ranges etc. will not be in- 
dicated again until required. The variational parameters 
x = {Di„} reside in the one-body density matrix p as 
given in Eq. (17). In order to work with familiar expres- 
sions in this section we will write p instead of x. 

Variation of the energy given in Eq. (18) with respect 
to p leads to the HF equations 

hnF[p] P = P h-n F [p] . (19) 
The one-body HF Hamiltonian 

hHF[p]=Y,\ l ) h ^[pU{k\ (20) 

i . k 

is a functional of the one-body density and its matrix 
elements are given by 



hj3F[p]ik 

or in short notation 

hwip] 



d 

Opk; 



£hf[c,/5] 
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(21) 



(22) 
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The eigenstatcs \4> v ) of /ihf [ P ] represent the basis in 
which both, /ihf[/5] and p are diagonal: 

hnF[p] = ^2 | <t>v ) £v ( 4>u | (23) 

p = J2\(f> v )n v ((f> v \ , (24) 

V 

where t v denotes the single-particle energy and n v are 
the single-particle occupation numbers that are zero or 
one in the case of a single Slater determinant. 

Eq. (19) represents the stationarity condition (5) 
and hence the HF approximation fulfills the Hcllmann- 
Feynman theorem. 

In nuclear structure theory the microscopic nucleon- 
nucleon interaction induces strong short-range correla- 
tions that cannot be represented by a single Slater de- 
terminant. Therefore the HF method as explained here 
cannot be used. For example the strong short-ranged re- 
pulsion makes all two-body matrix elements Vik.i m pos- 
itive and large, so that the HF Slater determinant does 
not give bound objects. The way out is to use effec- 
tive interactions that incorporate the short-range corre- 
lations explicitly, see for example Ref. [13-15]. Another 
approach is the density-matrix functional theory which 
we discuss in the following section. 



B. Hartree-Fock with density-matrix functionals 

It has turned out that bypassing the construction 
of an effective microscopic Hamiltonian by postulating 
an ansatz for the energy as functional of the one-body 
density-matrix p, as originally proposed by Skyrme for 
the non-relativistic and by J. Boguta and A.R. Bodmcr 
[16] and D. Walecka [17] for rclativistic nuclear physics 
(or by Kohn and Sham [18] for the atomic case), is very 
successful in describing ground state properties. The en- 
ergy functional £df[c,/S] contains a finite number of pa- 
rameters, c, that are adjusted by fitting observables to 
nuclear data. The shape of the functional is subject of 
past and present research [19-21] and is being improved 
to also apply for nuclei far off stability. Different from the 
HF case with Hamiltonian, Eq. (18), the densities may 
appear also with non-integer powers and the exchange 
terms are not calculated explicitly but absorbed in the 
form of the energy functional. 

Not all of the information residing in the one-body 
density-matrix p is used. Usually one uses the proton 
and neutron density p P (r), p n (r ), kinetic energy densities 
T p {r),T n (r), current densities j(f), etc. 

£df[c,p] = £ D f (c, p p (r ), p n (r),T p (r),T n (r), j (r), . . .) 

(25) 

In order to keep densities and currents consistent and 
corresponding to fermions they are expressed in terms of 



the single-particle states | <j> v ) 

1 0, > = £ 10 D ^ ( 26 ) 

i 

that represent the occupied states of a single Slater deter- 
minant. As they are expanded in terms of a working basis 
| i } the energy £df[c, p] is, like in the HF case, a function 
of the variational parameters x = {-D;„; i, v = 1, 2, . . . } 
or x = p. 

The difference to HF with a Hamiltonian is that the 
mean- field Hamiltonian /imf[p] obtained by the func- 
tional derivative of the density-matrix functional £ df 

/imf[/3] = tt£df[c,p] (27) 
bp 

is not given by a microscopic Hamiltonian anymore, but 
by the functional form of £df and the fitted parameters 
in the set c. 

The stationarity conditions (5) lead to the self- 
consistent mean-field equations 

/imf[/5] p = p Hmf[p\ (28) 

that have the same structure as the HF equations. 

Because the self-consistent solution is obtained by 
searching for solutions of the stationarity conditions (5). 
the Hcllmann Fcynman theorem (8) is fulfilled, even if 
one cannot refer to a microscopic Hamiltonian and a 
many-body state anymore. 

One should note that it is not mandatory that the 
single-particle states | (f> v ) with lowest single-particle en- 
ergies are occupied. Any combination of occupied states 
leads to a stationary solution fulfilling Eq. (28). 

Another interesting case is a one-body density with 
fractional occupation numbers < n v < 1 that may also 
commute with the mean-field Hamiltonian and hence ful- 
fills the stationarity condition (28) so that the Hcllmann- 
Feynman theorem is applicable. In such a situation one 
cannot attribute a single Slater determinant to the one- 
body density anymore, because p ^ p 2 is not idempotent. 

In any case the occupation numbers play the role of 
external parameters and should be regarded as members 
of the set c and not as variational parameters. 



C. Hartree-Fock-Bogoliubov 

The solution of the eigenvalue problem of /imf [ P ] pro- 
vides not only the coefficients of the occupied single- 
particle states but also a representation for empty states 
so that one has a complete representation of the one- 
body Hilbert space. With that one can define creation 
operators for fermions for occupied and empty states 

4 = J24 Diu , (29) 
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with 



,t 



) and |*)=cj|0> 



(30) 



and their corresponding annihilation operators a„ and Cj. 

Pairing correlations in the many-body state can be in- 
corporated by Bogoliubov quasi-particles that are created 
by 



al = u„al — v v a,v 
a% = u v a\ + v v a v 



(31) 



as linear combinations of the creation and annihilation 
operators, a v , dp, of the eigenstates of the one-body den- 
sity matrix (the so-called canonical states). The param- 
eters u v and v v can be chosen real and the requirement 
that a v and a\ are fermionic quasi-particle operators im- 



plies uf, - 



1. The pairing partner states v and v are 



usually mutually time-reversed states. 

The many-body trial state is expressed as 



*HF B )=n a i n ty i - ^ + 



0) , (32) 



where the product over fx runs over the so called blocked 
states or unpaired states and v runs over all other paired 
states. Besides the variational parameters residing in 
the operators aj, that create eigenstates of the mean- 
field Hamiltonian the energy depends now also on the 
variational parameters v v , hence x = u„; i, v = 

1,2,...}. 

As the trial state (32) has no sharp particle number 
the stationarity conditions Eq. (5) have to be augmented 
by a constraint on mean proton number Z and mean 
neutron number AT to obtain the selfconsistent Hartree- 
Fock-Bogoliubov (HFB) equations: 



£hfb = £ — \2, — \ n M 



(33) 



The proton and neutron chemical potentials, A p and A n , 
have to be regarded as members of the set c of external 
parameters. The additional constraints do not alter the 
arguments leading the Hcllmann Feynman theorem, thus 
it is also valid in the HFB case. 

In HFB it is convenient to introduce a generalized den- 
sity matrix 
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where p is the normal one-body density 
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and n the so called abnormal density 
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(34) 



(35) 



(36) 



Both can be expressed in terms of Di V and v v . The gen- 
eralized density matrix is idempotent and hermitian 



n 2 



n 



K and V) 
and the stationarity condition leads to 

h U f[k] ii = n Hmf[K] 



(37) 



(38) 



quite in analogy to the mean-field equations (28) without 
pairing correlations. The pseudo-Hamiltonian 



Hmf[^] = -~ ?HFB [c, H] 

51Z 



ft-MF — A 
—A* 
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MF 



(39) 



contains the mean-field Hamiltonian and the pairing part 
A. The chemical potentials A = (A p ,A Tl ) determine the 
mean proton and neutron number. For details and fur- 
ther reading see Refs. [22-24]. 

Let us write down only the one-body density because 
it is referred to in the results for 229 Th. The one-body 
density matrix is given by 



E 

v paired 



1+ E 

fj, blocked 



)n li (4> lt \, (40) 



where the | 0„ ) denote canonical basis states. The occu- 
pation numbers arc given by n v = rip = v 2 for the paired 
states and are the same for v and v. For the unpaired 
or blocked states n« = 0, 1. In the application we will 
consider only one blocked neutron state for 229 Th. 



IV. AMPLIFICATION 

To be more specific let us consider as external param- 
eters the fine structure constant a and the proton and 
neutron mass m p , m n , thus c = {a, m p , m„}. In this Sec- 
tion, we write down expressions pertaining to a Hamil- 
tonian; those characteristic to a density functional are 
entirely analogous, cf. Sec. II. The partial derivatives of 
a non-rclativistic Hamiltonian arc 



d 1 
—H(a,m p ,m n ) = — V c 
da a 

d „, 



dm r , 



z - 
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N - 
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m n 





(41) 
(42) 

(43) 



where Vc denotes the operator for the Coulomb energy, 
Z, N the proton and neutron number, respectively, and 
T p ,T n stand for the proton and neutron kinetic energy 
operator, respectively. 

According to the Hcllmann-Fcynman theorem a small 
variation 5a of the fine structure constant results in a 
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variation of an energy eigenvalue given by 



SE n 



dm r . 



(* n \V c \* n )+Za— ^ + Na 
da 
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( *„ | T n | „ ) 
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(44) 



In principle there could also be a dependence of the nu- 
clear interaction Vat on a, e.q., through meson masses, 
which we neglect here. The dependence of the nucleon 
masses on a can be estimated from the paper by Meifiner 
et al. [25] (Eq. (36)). Their estimate of the neutron- 
proton mass difference due to the electromagnetic inter- 



action is Al 



(EM) 



dm v 



-0.68 MeV which yields 



Am 



(EM) 
np 
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-0.36-10" 
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(EM) 
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(45) 



(46) 



A possible variation of some QCD constant c would 
lead to 

dm r 



dc 



N 



6E n =((* n \-^V N (c)\iS? n ) + Z 

dnip 

-(*n|T„|tf n ) 



-<*n|T p 



dm n 
dc 

dc 



(47) 



8c 



where Vn is the nuclear part of the interaction. This 
variation is actually linked to the dimensionless ratio c = 
m q / A-qcd, where m q denotes the current quark mass and 
Aqcd the strong interaction scale [8, 26, 27]. 

While the dependence of the nucleon mass on the cur- 
rent quark mass has been calculated [28, 29], the QCD 
constants enter the effective interaction Vn(c) in a very 
complicated and yet unknown way. In this paper we do 
not consider such variations of QCD constants explic- 
itly but calculate besides the total energies the kinetic 
energies and the Coulomb energies which then can be 
combined according to Eqs. (44) or (47) to obtain the 
variations with respect to variations of a or QCD param- 
eters. 

As the temporal variation of the fundamental con- 
stants a, nip and m n are tiny, if they exist at all, it has 
been proposed to consider transition frequencies 

w = E 1 (a) -E (a) (48) 

that can be measured with high precision. 

If the two energy levels belong to the same nucleus the 
terms with the proton and neutron number drop out and 
we obtain for the relative variation 8u> /cu of the transition 
frequency 

Slo 1 fdE 1 dE Q 



dc 



dc 



- AVb-AT, 
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da 



(49) 
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with the abbreviation 



AX 
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(50) 



for the difference of the expectation values of the oper- 
ators X = {Vc,T p ,T n } calculated with the stationary 
states. 

The results discussed in Sec. V show that AT n and AT p 
are of the same order as AVc so that the terms with the 
proton and neutron mass variations can be neglected. 

If there is a system, where the energy difference uj is 
much smaller than the difference between the Coulomb 
energies of the two states one gets an amplification factor 
in the measurement of 8a/ a. Flambaum [8] has proposed 
the nucleus 229 Th as a good candidate [7] because it pos- 
sesses two almost degenerate eigenstates with w w 8 cV 
according to a recent measurement [9] . 

The theoretical task is to investigate these states care- 
fully in order to get a reliable estimate for their Coulomb 
and kinetic energies. For the calculation of these quan- 
tities we employ in the following section state-of-the-art 
mean-field models and also include the effects of pairing 
correlations. 



V. RESULTS FOR 229 Th 

The nucleus 229 Th with 90 protons and 139 neutrons 
occurs in nature as the daughter of the a-decaying 233 U 
and decays itself with a half life of 7880 years, again by 
a emission. This nucleus has attracted lot of interest as 
it has the lowest lying excited state known. In Fig. 1 the 
spectrum is arranged in terms of rotational bands [30]. 
Two low lying rotational bands with = 5/2 + and 
3/2 + can be identified, with band heads that according 
to recent measurements differ in energy by only about 
8 cV. The first negative parity band with = 5/2 ~ 
occurs at 146.36 keV. 

Because of the amplification effect discussed in Sec. IV 
the transition 3/2 + — > 5/2 + is regarded as a possible can- 
didate to measure the time variation of the fine structure 
constant. A simple estimate of the moment of inertia for 
the two bands assuming a J (J + 1) energy dependence 
shows that the two rotational bands have very similar 
intrinsic deformation. Therefore, one does not expect 
a large difference between the Coulomb energies of the 
5/2 + and 3/2 + band heads that would be due to differ- 
ences in shapes. Instead, one has to consider and work 
out detailed effects due to different configurations of these 
states. 



A. Hartree-Fock with density-matrix functional 
theory 

To come to more quantitative statements we calculate 
the energies of these two states with the best methods 
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TABLE I: Total, Coulomb, neutron and proton kinetic energies of the 229 Th 5/2 + ground state calculated with different energy 
functionals. Differences of these energies between 3/2 + first excited state and 5/2 + ground state. 





Exp. 


SkM* 


SIII 


NL3 


5/2+ 


Ref. [30] 


HF 


HFB 


HF 


HFB 


RH 


E tot [MeV] 


-1748.334 


-1739.454 


-1747.546 


-1741.885 


-1748.016 


-1745.775 


V c [MeV] 




923.927 


924.854 


912.204 


912.216 


948.203 


Tn [MeV] 




2785.404 


2800.225 


2783.593 


2794.909 


2059.640 


T p [MeV] 




1458.103 


1512.705 


1442.018 


1477.485 


1106.697 


3/2+ - 5/2+ 
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AV C [MeV] 
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AT n [MeV] 
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AT P [MeV] 
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-0.022 
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FIG. 1: Measured low lying states of Th with spin and 
parity assignments [30]. Energies are in keV. 



available: density-matrix functional theory without and 
with pairing, and a spherical relativistic mean-field the- 
ory [31]. In the non-relativistic case we use the com- 
puter program HFODD (v2.33j) [32, 33] and employ two 
successful energy functionals, SIII [34] and SkM* [35]. 
Because of large uncertainties related to polarization ef- 
fects due time-odd mean fields [36], in the present study 
these terms in the energy functionals are neglected. The 
standard Slater approximation is used to calculate the 
Coulomb exchange energies. The code works with a 
Cartesian harmonic oscillator (HO) eigenstates as work- 
ing basis | i ) and allows for triaxial and parity-breaking 
deformations, but the states determined in 229 Th turn 
out to have axial shapes with conserved parity. In our 
calculations we have used the basis of HO states up to 
the principal quantum number of Nq = 18 and the same 



HO frequency of Huj = 8.05 MeV in all three Cartesian 
directions. 

In the discussion, the resulting basis | </>„ ) will be as- 
signed Nilsson quantum numbers f2 7r [N, N z , A] (for de- 
tails see Ref. [22, 37]) by looking for the largest over- 
lap \(n 7r [N,N z ,A] 1^)1 with a Nilsson state. TV = 
N x + N y + N z = N p + N z + A denotes the total number of 
oscillator quanta, N z the number of quanta in the direc- 
tion of the symmetry axis. A = N p , N p — 2, ... or 1 and 
Q = A ± i are the absolute values the projection of or- 
bital angular momentum and total spin on the symmetry 
axis, respectively. 7r = (— 1) N is the parity. 

It turns out that for both functionals the lowest HF 
state has f2 7r = 5/2 ~ and thus corresponds to the intrin- 
sic state of the K n = 5/2~ band that is experimentally 
located at 146.36 keV. The total binding energy amounts 
to -1739.454 MeV for SkM* and -1741.885 MeV for SIII 
which should be compared to the experimental energy of 
-1748.334 MeV [30]. On this absolute scale the calcu- 
lated HF energies are already amazingly good. 

In order to get the experimentally observed parity we 
perform the variation procedure in the subspace with 
positive parity. For both energy functionals the Slater 
determinant with the lowest energy has = 5/2 + . All 
energies are summarized in Table I. The last neutron 
occupies the level labeled with 5/2+ [622] in the SkM* 
and with 5/2 + [633] in the SIII case. One should keep in 
mind that the single-particle states are superpositions of 
Nilsson states and the labeling refers only to the largest 
component. The states are detailed in Eqs. (51)- (54) be- 
low. 

The neutron single-particle energies e v and the occu- 
pation numbers arc displayed in Fig. 2(a) for SkM* and 
in Fig. 3(a) for SIII. The 5/2" [752] state which would be 
occupied in the negative parity case is very close to the 
5/2 + [633] state - for SIII they are almost degenerate. 

As can be seen from Table I, the total HF binding 
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FIG. 2: (Color online) Neutron single-particle energies and occupation numbers labeled by the asymptotic Nilsson quantum 
numbers v = [N , N z , A] for the SkM* energy functional, (a) HF mean- field energies e„ and occupation numbers n v . Full 
bars denote n v + — 2, i.e. two particles in degenerate pair of states with nij = ±fi. Half full gray (pink) bars denote one 
particle, n M = 1, np = 0. (b) HFB mean-field energies e„ (eigenvalues of /imf) and occupation numbers vt. Length of bars 
indicates n v + n^. Gray (pink) bars stand for blocked states with n M = 1, np = 0. 
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FIG. 3: (Color online) Same as Fig. 2, but for the SIII energy functional. 



energy agrees with the measured one up to about 9 MeV 
for the SkM* and up to about 6 MeV for the SIII energy 
functional. Keeping in mind that no parameters have 
been adjusted to the specific nucleus considered here it is 
surprising that these mean-field models can predict the 
energy with an uncertainty of only about 0.5 %. 

Putting two neutrons in the 5/2+ [622] (for SkM*) or 
5/2+ [633] (for SIII) level and one in the 3/2+ [642] (for 
SkM*) or 3/2+[631] (for SIII) level and minimizing the 
total energy yields an excited HF state that is to be re- 
garded as the intrinsic state of the experimentally ob- 
served = 3/2+ band. As can be seen in Table I, 



the excited states occur at 0.619 MeV for the SkM* and 
at 0.141 MeV for the SIII density functional. The dif- 
ference in Coulomb energies AVc amounts to 0.451 MeV 
for SkM* and to -0.098 MeV for SIII. The kinetic energy 
differences AT p and AT n for protons and neutrons, re- 
spectively, dissent even more. These deviations between 
the two energy functionals reflect the differences in the 
structure of the intrinsic states as also seen from the dif- 
ference in the single-particle states discussed above. 
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B. Relativistic mean- field 

We also perform a spherical relativistic mean-field cal- 
culation with the NL3 parameter set [38] and find that in 
the ground state the last neutron occupies a 2g g / 2 orbit. 
For vanishing deformation the Nilsson state 5/2 + [633] 
belongs to the subspace spanned by the spherical 2g 9 / 2 
orbits so that this result is not unreasonable when com- 
paring with the deformed SIII calculation. The leading 
component for the 3/2 + state is for SIII the 3/2 + [631] 
orbit (cf. Fig. 3) which for deformation zero belongs to 
the l«ii/2 subshell. Therefore we create the excited state 
with SI* = 3/2+ by a particle-hole excitation from the 
l*n/2 to the 2^9/2 shell so that there are 11 neutrons in 
the lin/2 and 2 neutrons in the 2g 9 / 2 shell. 

The resulting energies are listed in the last column 
of Table I. The total binding energy is similar to the 
non-relativistic one, but the particle-hole excited state 
is 2.41 MeV higher. Different from the deformed mean- 
field the single-particle states of a spherical potential have 
good total spin j and are (2j + l)-fold degenerated with 
large gaps between them. The single-particle energy dif- 
ference £9/2 — £11/2 = 2.74 MeV explains the large excita- 
tion energy of the particle-hole pair of 2.41 MeV, which 
includes the rearrangement energy. 

We conclude that a spherical calculation is not appro- 
priate for this particular question. In a recent publication 
[39] a similar spherical relativistic mean-field calculations 
comes to results comparable to ours with a Coulomb en- 
ergy difference of 0.7 MeV. Because of the unphysical 
properties of a spherical 229 Th we will not commit our- 
selves to the spherical case any longer but proceed to 
consider the effects of pairing. 

C. Hartree-Fock-Bogoliubov 

We include the pairing correlations with the Bogoli- 
ubov ansatz (32) and perform a self-consistent HFB cal- 
culation based on the SkM* and SIII density-matrix func- 
tional. For the SIII case proton and neutron pairing 
strengths of Vb = —260 and — 180 MeV fm 3 (for a volume- 
type contact force) are adjusted to reproduce the total 
binding energies and the odd-even staggering with the 
neighboring nuclei. Proton and neutron density ma- 
trices and pairing tensors are calculated by including 
contributions from quasiparticle states up to the cutoff 
energy of 60 MeV. Calculations are performed by self- 
consistcntly blocking the 5/2 + and the 3/2 + quasipar- 
ticle states. In the 5/2 + and 3/2 + configurations, this 
yields the average HFB proton and neutron pairing gaps 
of A p = 2.4 MeV and A„ = 0.65 MeV, and A p = 2.4 MeV 
and A„ = 0.68 MeV, respectively, for the SIII case and 
slightly larger values for the SkM* case. 

The results for the energies are summarized in Table I 
and for the mean-field single-particle energies and the 
occupation probabilities in Fig. 2(b) and Fig. 3(b). The 
first to note is that the excitation energy is improving. 



Its value decreases from 619 keV down to —46 keV for 
SkM* and from 141 keV down to -74 keV for SIII. On the 
accuracy level one can expect from this model the 5/2+ 
and the 3/2 + states are degenerate, like in experiment. 

By looking at the occupation numbers displayed in 
Figs. 2(b) and 3(b) one sees that about 5 single particle 
levels near the Fermi edge assume fractional occupation 
numbers significantly different than or 2. In the SIII 
case (Fig. 3(b)) they are almost identical for the 5/2+ and 
3/2+ states except for the 5/2+ [633] and 3/2+ [631] levels 
that switch their role. For the SkM* case (Fig. 2(b)) the 
blocked states are energetically further away from each 
other which causes more deviations in the occupation 
numbers. 

This characteristic pattern of occupation numbers ren- 
ders the HFB results qualitatively different than the HF 
ones. Indeed, in the HF case either the 5/2 + or the 3/2 + 
orbital has the occupation number equal to 1. Therefore, 
polarization effects exerted by these two orbitals are able 
to render different values of observables calculated for 
the 5/2 + and 3/2 + states. In the HFB case, occupation 
numbers of the 5/2 + and 3/2 + orbitals are close to 1 for 
both 5/2 + and 3/2 + configurations. Differences in ob- 
servables may here only occur due to the fact that the 
occupation number of the blocked state is exactly equal 
to 1 , while that of the other state is approximately equal 
to 1, depending on its closeness to the Fermi level. 

In the paired case (HFB), we are faced with the sit- 
uation, where the zero-order approximation renders ob- 
servables calculated in the 5/2 + and 3/2 + exactly equal. 
Indeed, such equality would be the case for the 5/2 + 
and 3/2 + orbitals located exactly at the Fermi surface 
and having exactly the same occupation numbers. Note 
that values of all observables calculated with the HFB 
approach depend only on the canonical states and oc- 
cupation numbers, irrespective of which state has been 
blocked in obtaining them. Of course, one can obtain 
different occupation numbers for both orbitals in ques- 
tion when they are energetically split. However, this may 
contradict the experimental energetic degeneracy of the 
corresponding configurations. All in all, within a zero- 
order paired approach, polarization effects of the 5/2+ 
and 3/2 + orbitals become exactly averaged out and the 
anticipated differences in observables can occur only due 
to first-order corrections. 

This fact is perfectly well visible in our results for the 
Coulomb energy differences shown in Table I. For the 
SIII energy functional, only 1 keV remains for A Vc ■ This 
reduces the amplification factor of Eq. (49) to about 100. 
For SkM* a larger value of AV C of about 300 keV is 
obtained due to a larger splitting of the corresponding 
single-particle orbitals. From this one must conclude that 
pairing correlations result in two states with even more 
similar charge distributions than in the HF calculation. 
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TABLE II: Rms-radius and intrinsic quadrupole moments of neutron and proton densities of the 229 Th 5/2 + ground state 
calculated with different energy functionals. Differences of these moments between 3/2 + first excited state and 5/2 + ground 
state. 





SkM* 


SIII 


c; /9+ 


1 1 1 


rir d 


1 1 1 


rlr id 


Rrms (neutron) [fm] 


5.8789 


5.8716 


5.8971 


5.8923 


Rrms (proton) [fm] 


5.7180 


5.7078 


5.7817 


5.7769 


Q20 (neutron) [fm ] 


9.4407 


9.2608 


9.1990 


9.0711 


Q20 (proton) [fm 2 ] 


9.5461 


9.3717 


9.3542 


9.1643 


3/2+ - 5/2+ 










A Rrms (neutron) [fm] 


-0.0040 


0.0036 


-0.0008 


-0.0005 


AR rms (proton) [fm] 


-0.0038 


0.0039 


0.0000 


-0.0005 


AQ20 (neutron) [fm 2 ] 


-0.2427 


0.2647 


-0.0767 


-0.0516 


AQ20 (proton) [fm 2 ] 


-0.1824 


0.2756 


-0.0339 


-0.0495 



D. Radii and quadrupole moments 



In Tabic II the rms-radii and quadrupole moments 
(normalized with proton and neutron number, respec- 
tively) of the neutron and proton point-densities are 
given for the ground state and as differences for the ex- 
cited state. The quadrupole moments of the protons are 
somewhat larger than for the neutrons. When reoccu- 
pying the last neutron the neutron quadrupole moment 
decreases substantially in the SkM*-HF calculation and 
drags along via the nuclear interaction the quadrupole 
moment of the protons. At the same time both rms- 
radii are decreased. A smaller charge radius and smaller 
charge quadrupole moment are consistent with the in- 
crease of the Coulomb energy which explains the large 
positive AVc in the SkM*-HF case. When including 
pairing the effect goes in the opposite direction for SkM*- 
HFB. 

For the SIII functional HF and HFB calculations do 
not lead to noteworthy changes in the moments when 
reoccupying the last neutron and thus the Coulomb dif- 
ferences remain also small. This can be anticipated when 
looking at the single-particle energies of the involved neu- 
tron states and the occupation numbers (Fig. 3). 

The mean-field single-particle state that corresponds 
to the blocked HFB state occupied by the unpaired neu- 
tron is represented in Nilsson orbits for SkM*-HF as 



and for SkM*-HFB including pairing as 
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3 + 



.509| 622 > 
.402 j 633) 

.010| 622) 
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.266(862) 
.249(862) 



(51) 



+ > 



+ .504(622) +.487 
+.418(633) -.383 
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The SIII-HF calculation gives 
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and the SIII-HFB with pairing 
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642) 
613) 

642) 
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-.159(862) 
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In SIII-HF and SIII-HFB the blocked states exhibit more 
concentration on the dominant Nilsson orbits [633] and 
[631]. As both have the same nodal structure in z- 
direction (N z = 3) one expects less difference in the den- 
sity distribution than in the SkM* case where both states 
are superpositions of several Nilsson orbits with similar 
amplitudes. In the SkM* case the leading component of 
I I + ) is [622] and of | | + ) is [642], which implies a change 
in nodal structure in z-direction. In summary the polar- 
ization of the proton distribution due to the reoccupation 
of the level with the unpaired neutron has less effect in 
the SIII than in the SkM* case. 

In Ref. [40] the finite-range microscopic-macroscopic 
model has been used to study the problem. The authors 
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find small Coulomb energy differences similar to our SIII 
case. Also the decomposition of the last neutron orbits 
[47] shows a mixture of several Nilsson orbits resembling 
more our SIII states than the SkM* states. 

We should like to point out that the Nilsson orbit 
5/2 + [633] that is usually used to classify the = 5/2 + 
ground state band [9, 41, 42] does not even contain a 
single-particle spin j 7 ' = 5/2 + . Due to positive parity 
and A = \mi\ = 3 the orbital angular momenta con- 
tained in 5/2 + [633] arc I = 4,6,.... Thus the lowest 
possible spin in the 5/2 + [633] Nilsson state is j 7r = 7/2 + . 
This implies that in a core plus valence-neutron picture 
the 5/2+ [633] state needs to be coupled to the excited 
jit _ 2+ state in 228 Th in order to get the ground state 
spin J* = 5/2 + . In the deformed mean- field description 
the total angular momentum J T = 5/2+ of the nucleus 
arises from both the 5/2 + [633] orbital and the underlying 
deformed 228 Th core. 

In Rcf. [41. 43] experimental data have been 
compared with structure calculations for 229 Th that 
use the quasiparticle-phonon model [44] employing a 
phcnomcnological Nilsson mean field and multipole- 
multipole residual interactions. A very good description 
of the 229 Th level structure has been achieved by an ap- 
propriate fit of the interaction parameters. It has been 
found that the coupling of the single-quasiparticle de- 
grees of freedom to the collective octupole vibrational 
state of the 228 Th core is essential to reproduce the par- 
ity partner bands observed in experiment. However, as 
this model is not self-consistent it can not be used for our 
considerations. 

In a fully consistent calculation scheme based on a rcl- 
ativistic density-matrix functional [45] it has been shown 
that the coupling to low lying vibrations noticeably im- 
proves the description of the single-particle spectrum 
around the Fermi surface including the ordering of the 
levels. This model is however programmed only for spher- 
ical cases. 

One has to realize that the situation is not so sim- 
ple, valence and core nucleons have to be considered self- 
consistcntly and in the next generation of nuclear models, 
besides the projection of the deformed intrinsic state on 
total spin and particle number, one also may need to go 
beyond the mean-field picture by coupling to low-lying 
core excitations. 



E. Predictive power and accuracy of observables 

In the energy-density functional picture one gives up 
the explicit knowledge of a microscopic Hamilton opera- 
tor H(c) acting in many-body Hilbert space. Its expec- 
tation value is replaced by the energy functional £ [c, p] 
from which one cannot refer back to the Hamiltonian. 
It is important to note that one can also not refer back 
to a many-body state | ^ } that represents an approx- 
imation to a true stationary eigenstate of H(c). One 
can of course construct a single Slater determinant with 



the operators aj, which create the occupied states, but 
this Slater determinant is more an auxiliary object that 
ensures quantum properties like Pauli principle or uncer- 
tainty relation. This Slater determinant misses for ex- 
ample various kinds of typical nuclear correlations that 
exist in the true eigenstate. 

This raises the important question if one has predic- 
tive power for other observables than the energy. The be- 
lieve is that observables which can be calculated from the 
one-body densities that appear in the energy functionals 
should be trusted. In our case we calculate Coulomb and 
kinetic energies, which are given by densities that are 
included in the set of variational variables of the energy 
functional £[c, p] and therefore should be predicted with 
high accuracy. 

Besides these more general considerations there are 
also concerns about numerical precision. In the SkM* 
case the Coulomb energy difference AVc = —0.307 MeV 
= (924.854-925.161) MeV is a result of subtracting two 
big numbers that have been calculated numerically. That 
means a precision of 10 keV for each of them is desirable. 
For SIII-HFB, where AVb=0.001 MeV, a precision of 
0.1 keV is needed. We have checked that we can reach 
enough numerical precision by sufficient iteration steps. 

But there is also the quest for accuracy. Let us for 
example consider the approximate treatment of the ex- 
change term. Its contribution in HFB-SIII is —34 MeV, a 
10% error means already 3 MeV uncertainty in ( Vc ) = 
y direct _^yexch _ Q n ^ e other hand one would expect that 
this is mainly a systematic error that is similar for the two 
states so that the difference should be less affected, but 
1% still means 0.3 MeV error. The actual situation for 
SIII-HFB is AVg irect = 0.29 keV and AV§ xch = 0.71 keV 
which adds up to the AVb=0.001 MeV listed in Ta- 
ble I. In this case the approximate exchange term gives 
the larger contribution which weakens strongly the con- 
fidence in the calculated value of the amplification. 



TABLE III: Violation of the Hellmann-Feynman theorem S = 
a (dE tot (a )/da\ -E Coul (a a ) in HF and HFB calculations 
with SIII Skyrme functional for 228 Th. 



HF 


HFB (no cutoff) 


HFB (with cutoff) 


0.02 keV 


0.5 keV 


200 keV 



As discussed in Sec. II the proper derivation of the 
equations of motion from an energy functional is cru- 
cial for the validity of the Hcllmann Fcynman theorem. 
Without this theorem one would calculate numerically 
the derivative of the energy with respect to a and cre- 
ate a new source of numerical errors. We tested the va- 
lidity of the Hellmann-Feynman theorem by comparing 
to numerical derivatives using a five-point formula. In 
Table III the deviations are listed for the 228 Th ground 
state of the SIII functional. In the HF calculation the 
deviation is within the numerical uncertainty induced by 
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the five-point formula. In the HFB calculation one gets 
also sufficient accuracy when no cutoff in the quasipar- 
ticle subspace contributing to the pairing interaction is 
applied. But the accuracy drops by three orders of mag- 
nitude when the density matrices and pairing tensors are 
calculated by including contributions from quasiparticle 
states up to the cutoff energy of 60 MeV. The reason is 
that this truncation does violate the variational structure 
of the HFB equations [46]. For most observables, this vi- 
olation induces small effects, and usually can be safely 
neglected, but it does show up in the very demanding 
calculation of the Coulomb energy differences. 

As discussed before, another accuracy issue is the fact 
that our model does not contain projection on good total 
spin and sharp particle number. Also the possibility that 
configurations could admix that consist of the K v = 1 ~ 
band of 228 Th coupled with a single-neutron 5/2~ state 
cannot be excluded. These questions have to be the task 
of future investigations. 

VI. SUMMARY 

We have investigated the lowest two states of 229 Th 
that are almost degenerate in energy. Two very success- 
ful energy functionals SkM* and SIII have been employed 
in a density-matrix functional theory. Hartree-Fock 
and Hartree-Fock-Bogoliubov calculations have been per- 
formed and compared. The result is that for the SkM* 
functional the difference in Coulomb energy AVc be- 
tween excited and ground state ranges from 450 keV 
without pairing to —300 keV when pairing effects are 
included. On the other hand the SIII-HFB result gives 
AVc — 1 keV only. The differences in neutron and pro- 
ton kinetic energies are of similar size and also quite dif- 
ferent for SkM* and SIII. 

Altogether, the nuclear models we used predict am- 
plification factors A = AVc/u>, between the drift of the 
transition frequency Suj/lo and the drift 8a/ a in the fine 
structure constant, that have absolute values varying be- 
tween about 10 2 and 10 4 . We have pointed out and 
discussed the fact that the pairing correlations smooth 
out polarization effects exerted by the single-particle or- 
bitals. Therefore, such correlations not only dramatically 
decrease the anticipated amplification factors but also 
make their determination very uncertain, due to depen- 
dence on very detailed properties of the mean-field and 



pairing effects. 

We have also performed spherical calculations and con- 
clude that spherical models should not be consulted as 
they are too far from reality to provide serious numbers. 

As even the sign of the amplification factor is uncer- 
tain, much more refined calculations arc needed that 
include coupling to low-lying core excitations and pro- 
jection on eigenstates with good total angular momen- 
tum and particle number. Before being able to provide 
reasonably trustable numbers how the transition energy 
varies as function of the fine structure constant a one 
has to make sure that the model reproduces the three 
low lying rotational K 1 * = 5/2 + , 3/2 + , 5/2~ bands up 
to J 9/2 and the known electromagnetic transitions 
within the bands and between them. This would provide 
more confidence in the quality of the many-body states 
and their Coulomb energy. 

In any case the calculations must treat all nucleons (no 
inert core) because the whole effect comes from a subtle 
polarization of the core protons. Furthermore the model 
has to be of variational type in order to make use of the 
Hellmann-Feynman theorem. Without that one cannot 
be sure that the polarization effects caused by the strong 
interaction are treated consistently with the necessary 
accuracy. Rough estimates and simple minded models 
are not sufficient. 

The experimental endeavor for measuring the drift of 
the transition frequency in 229 Th has to be accompanied 
by substantially improved models on the nuclear theory 
side and attempts to gain experimental information on 
the Coulomb energies via radii, quadrupolc moments or 
even form factors. Without a concerted action of experi- 
mental and theoretical efforts the goal of improved limits 
on the temporal drift of fundamental constants cannot 
be reached. 
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